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§^' 1 Introduction 

^ . QCD, in contrast to QED, is not exhausted by perturbation theory. While 

QCD perturbation theory appears to be very successful at higher momentum 
>■ I transfers where the strong coupling constant becomes small, some of the 

most fundamental properties of hadronic physics like color confinement and 






O ' hadronization of quarks and gluons, characteristic for the low momentum 

^ . transfer region, are not accessible to perturbative methods. 



O 



O ■ QCD is a nonabelian gauge theory. In 1975, Polyakov et al. |]T| made the im- 

portant discovery that such theories in their pure gauge form, i.e. without the 
involvement of spontaneous symmetry breaking, have localized classical solu- 
tions in Euclidean space-time. Since the objects described by these solutions 
are localized in 3-1-1 dimensions, they describe "events" rather than normal 

>-Ch ■ particles which are localized in 3 dimensions. Accordingly, they were named 

pseudoparticles by Polyakov [1] and instantons by 't Hooft §. The replace- 

K> ! ment of time t by ix^ (real X4) by going into Euclidean space corresponds 

H I to the method of calculating tunneling amplitudes in the nonperturbative 

WKB approximation, where the characteristic exponential suppression fac- 
tor exp(— 87r^/(7^) is obtained from solving the classical equations for imag- 
inary time. Instantons are, therefore, interpreted as tunneling events. The 
tunneling, of course, takes place in real (Minkowskian) space, although it is 
actually computed in Euclidean space. As was first realized by Ringwald and 
Espinosa 0, the exponential suppression of tunneling processes may be over- 
come at high energies by multiple emission of gauge bosons accompanying 
fermion production. 

The tunneling events violate a conservation law that is obeyed by perturba- 
tive solutions: For the classical (BPST) solutions there exists a charge like 
quantity 
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S27r2 
d^xF,,F,, = ±^ (1) 



where F^^ is the gauge field strength tensor and F^y its dual tensor. Since 
from the Adler-Bell-Jackiw axial triangle anomaly 0] {g"^ /lQ'K'^)F^yF^y equals 
the anomalous divergence of the axial- vector divergence, tunneling events i.e. 
instantons connect states of different axial charge Q^if) = J d^xJ^{x,t) and 
one obtains 

/ + 00 P 

dt / d^xdoJoix) = ±2. (2) 

-oo J 

This relation holds separately for each fiavour and for Nf fiavours one has 

AQs = ±2Nf. (3) 

Due to their topolocical nature (1) and (2) hold in Euclidean as well as in 
Minkowskian space-time. The violation of the conservation of axial charge 
i.e. helicity is a most significant property of instanton events. 

The relevance of instantons for different aspects of hadronic physics has been 
reviewed recently |^. In the following, we discuss recent work on the role of 
QCD-instantons in deep inelastic electron proton scattering. 

2 Instantons in Deep Inelastic Electron Pro- 
ton Scattering 

Balitsky and Braun ||^ first showed that the contribution of instanton (/)- 
induced processes to deep inelastic electron proton scattering from a real 
gluon calculated in Euclidean space and continued to Minkowski space, rises 
very rapidly with decreasing Bjorken-a;. Their calculations were restricted, 
however, to a; > 0.3 — 0.35. A detailed theoretical |]^, ||, ^ and phenomeno- 
logical [jlO|, |Tl| investigation of deep inelastic scattering at HERA has been 
pursued by Ringwald, Schrempp and collaborators. 

The leading /-induced process contibuting to deep-inelastic ep-scattering, 
given by a hard quark originating from the virtual photon and reacting in 



the presence of an instanton with a gluon from the proton, is shown in fig. 
1. 




Figure 1: Leading graph for instanton- 
induced contribution to deep inelastic 
electron proton scattering. It denotes 
the instanton subprocess with the ob- 
servables Q'"^ = —q''^ = —{q — q')'^, 
x' = Q'^/2{g . q'), Wf = {q' + gf = 
Q'^{l-x')/x'. 



The corresponding contribution with a quark originating from the proton is 
suppressed by a1 and can be neglected in the kinematical region of interest 
(see below), where the gluon density is high. 

The inclusive deep-inelastic ep cross section in /-perturbation theory in the 
Bjorken limit can be approximated as 
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(4) 



where dCgg/dx'dQ''^ is a differential luminosity corresponding to a convolu- 
tion-like integral over photon flux, quark flux and gluon density and cr^ is 
the total cross section of the /-induced quark-gluon subprocess. 

The critical part of (4) is the calculation of aj, which contains all /-dynamics. 
Outlining the essential dependencies, it can be written as [B], 01 
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where the integration is performed over the "collective coordinates" of the 
instanton i.e. the size p weighted with the density distribution function D{p), 
the / — /-distance vector R and the /- /-relative colour orientation U. I as 
well as / enter here since the cross section results from the modulus squared 
of the amplitude for processes induced by single / resp. /. The function 
Q with the large factor —An/as contains multiple emission of gluons. D[p) 



is known in /-perturbation theory for a{fir) ln(p/ir) ^ 1 where /i^ is the 
renormahzation scale. D{p) has a powerlaw behaviour. 

D{p) ~ p6-|„,+0{a) (g) 

which in general results in a breakdown of /-perturbation theory for large 
values of p. In deep inelastic scattering, however, the exponential term [7] 
containing Q' supresses the contribution form large size instantons and is the 
key for the applicability of /-perturbation theory. 

For calculating the cross section for deep inelastic scattering DIS at HERA, 
the range of collective coordinates has to be known in which /-perturbation 
theory can be relied on. This range of collective coordinates then has to be 
translated into the relevant DIS-variables. 

For estimating the range of validity of /-perturbation theory QCD lattice re- 
sults from the UKQCD collaboration [12] have been extremely useful. Fig. 2a 
shows a comparison of /-perturbation theory with lattice-QCD in the con- 
tinuum limit [9] for the /-size distribution. The normalization is practically 
independent on the renormahzation scale Pr, but strongly dependent on the 
value of Ajfg, for which the recent accurate result from the ALPHA collab- 



oration |T^ ^M5(nf=o) ~ (238 ± 19)MeV has been taken. It has to be noted 



that the pertubative result is essentially parameter free. Fig. 2b shows the 
corresponding comparison for the / — /-distance contribution. Very good 
agreement down to / — /-distances R/ < p~ >^ 1 is observed. 

In the deep inelastic regime the integrals over collective coordinates are dom- 
inated by a unique saddle point p| which relates the collective coordinates 
to the Bjorken variables Q', x' of the /-subprocess as 

p\-p*^ 1/Q', R*^ ~ l/{p + q'f ^ R*^/p*-p* ~ g'V(p + q'f = x'/l - x'. 

(7) 
The limits on the collective observables then transform into limits of the DIS 
variables in the following way: 

p* < 0.30 - 0.35/m Q'/^ws > 30-8 

(8) 
R*/p* > 1 x'> 0.35 

With these cuts, additional general cuts of x > 10~^ and 0.1 < y < 0.9 
and a "technical cut" of Q'^m = Qmin^ i^ order to suppress contributions 
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Figure 2: Comparison of "collective coordinates" obtained from instanton- 
pertm'bation theory (cm've) with the continuum limit of equivalent lattice QCD 
data (points): a) /(I)-size distribution; b) / — /-distance distribution, normalized 
to the value at R= 1.5/m. From Refs. g, ^ 

from nonplanar graphs which are hard to calculate, a total cross section 
for /-induced events of 29 ± 10 pb is obtained, where the errors refer to the 
uncertainty in Aj^^ only. For a total integrated luminosity of about 100 pb~^, 
which will be reached for each of the HERA experiments HI and ZEUS by fall 
2000, this corresponds to about 3.000 /-events. Since the total cross section 
for standard DIS events is several orders of magnitude higher it is evident, 
however, that a search of /-induced effects has to make use of specific event 
characteristics. 



Search Strategies for /-induced Events at 
HERA 



The Monte Carlo generator QCDINS |T^, developed on the basis of the 
perturbative model discussed above for the detailed simulation of /-induced 



events, yields the following event characteristics |TU[: 



a jet ("current jet") from the outgoing quark in the primary photon 
interaction (compare fig. 1) 

after removal of the current jet a uniformly populated band of hadrons 
in rj — 0-space (with the pseudorapidity rj and the azimuthal angle (p) 
resulting nearly in isotropic decay in the hadronic CMS (defined by 
q + P = 0) 

high parton multiplicity (for 3 flavours: 3 (gg)-pairs; at a^ a mean value 
oiO{l/as) — 3 gluons) resulting in a mean charged hadron multiplicity 
of about 20 

"flavour democracy" resulting in a significantly higher fraction of strange 
hadrons in /-induced than in standard DIS-events 

a higher total transverse momentum of ~5 GeV/rj-nnit as compared to 
~ 2 GeV/?7-unit for DIS-events |l^ 



A typical MC-event is shown in fig. 3. 
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Figure 3: Distribution of the trans- 
verse energy Et in pseudorapidity 
(r/)-azimuthal ((/))-plane in the hadronic 
CMS for a typical instanton-induced 
HERA-event generated by QCDINS 
(x = 0.0012,Q2 = 66 GeV2,pT(Jet) = 
3.6 GeV) after typical detector cuts. 
Clearly recognizable are the current jet at 
(j) = 160°, ry ~ 3 and the instanton band 



at ~ 7? ~ 2. From Ref. Q 



Search strategies for /-induced events in a background of normal DIS events 
based on Monte Carlo data have been discussed by Carli, Gerigk, Ringwald 
and Schrempp [0. The aim is to enrich /-induced events in a data sample 



by cutting on selected observables while optimizing the separation power, 
defined as the ratio of the detection efficiencies for /-induced and DIS-events; 
in addition, a lower limit of 10% for the /-detection efficiency is required. 
Fig. 4 shows the optimization for the number of charged hadrons contained 
in the /-band after removal of the current jet. 
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Figure 4: a) Normalized MC-generated distributions of the number of hadrons n^ 
in the instanton band for instanton induced (INS) and normal events (DIS) after re- 
moval of the current jet hadrons. b) Maximum separation power INSeff/DISeff 
as function of cuts on n'^ (see (a)) keeping the instanton efficiency INS^jf at 



> 10%. From Ref. |ll| 



Cuts on single observables typically result in a separation power of ~ 20 and, 
therefore, cuts on several observables have to be combined in order to achieve 
a higher background reduction. Correlations among the different variables, 
however, tend to weaken the effect of combinations. For a set of 6 selected 
observables, which are shown in fig. 5, a separation power of about 130 was 
obtained. 



4 Recent HI Results 



Based on this experience, the technique for enhancing /-induced events has 



recently been applied by the HI collaboration to real data ||T6| : Out of a set 
of 6 observables a subset of 3 observables has been selected on which cuts 
are applied (called in the following "primary observables") while the effect 
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Figure 5: Normalized distributions of characteristics observables for normal DIS 
events (ARIADNE) and /-induced events (QCDINS + HERWIG). Shaded bands 
indicate model uncertainties, arrows show optimized cuts with the allowed region 
left of the cut. Variables: pT and Q'^ refer to the current jet (compare fig. 1), 
E'rp ^ and n'^ are the total transverse energy and number of hadrons, resp. in the 
/-band after removal of the current jet, A is a measure of the isotropy of events 



(see text), //lo are the Fox- Wolfram-moments. From Ref. |11] 



of these cuts on the distributions of the other observables not subjected to 
cuts (called "secondary observables") is studied. The preliminary results 
are based on a total integrated luminosity of 15.8 pb~^ collected by HI in 
1997 and corresponding to about 260 k DIS-events within fixed cuts of x > 



8 



10-3, 0.1 < y < 0.6 and Q^ > 156°. 

The primary observables chosen for the analysis were the number of hadrons 
rih in the J-band, the sphericity Sph of the event in the hadronic CMS (both 
after removal of the current jet) and the squared momentum transfer Q'^ of 
the /-subprocess (c.f. fig. 1). The normalized distribution of these observ- 
ables are compared in fig. 6 with two quite different models for DIS processes, 
i.e. the colour dipole model CDM |l^ and the matrix element-parton shower- 
model MEPS |8[ as implemented in the Monte Carlo ARIADNE resp. RAP- 
GAP. Agreement on the level of 5-20% is observed. The signal expected for 
/-induced events from QCDINS is 10^ — 10^ below the DIS background. 

A systematic investigation of the effect of cutting primary observables has 
been carried out. The following combinations of cuts have been studied: 
n6>5,6,7,8,9; Sph > 0.4,0.5,0.55,0.6,0.65; Q''^{GeV^): 95,100,105,110, 
115 < Q'^ < 200. The results have been classified in 3 scenarios: "A" for the 
highest instanton efficiency eiNs, "B" for high instanton efficiency at good 
separation power eiNs/^Dis and "C" for the highest separation power at an 
instanton efficiency > 10% (see table 1). 

Table 1: 
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95.0-200.0 


0.40 
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32% 
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0.40 
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21% 
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105.0-200.0 


0.50 
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11% 


86 
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As example, the effect of cut C on primary observables is shown in fig. 7. 
A background reduction by a factor of 600-800 is achieved. The remaining 
549 events are to be compared with 362 dz 25 events expected from CDM 
and 435±30 events from MEPS. For the resulting distributions the shape of 
the excess with respect to the DIS models is compatible with the QCDINS 
expectation. It has to be noted, however, that the excess is of the same 
magnitude as the difference between the two DIS models and that also the 



QCDINS prediction is subject to systematic uncertainties [|T9|. 

As secondary observables were chosen the total transverse energy Efb of 
hadrons contained in the J-band, the transverse energy Etjet of the current 
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Figure 6: Observables used to cut ("primary observables" ) : Comparison of data 
with Monte Carlo (see text). From Ref. [16] 



jet and a third quantity A, 

m 



m 



which is a measure for the isotropy of events 
It is obtained by finding the two axes with respect to which the pro- 
jection of all hadronic momenta (after removal of the current jet) is minimal 
(Ein) resp. maximal (Eout) and forming {Ein — Eout) / Ein = A^. A;, measures 
the i?rweighted azimuthal isotropy; it takes small values for isotropic events, 
expected for /-induced reactions, and large values for pencil-like 2-jet-events, 
expected to constitute a major fraction of the DIS-background. 
The effect of cut C on the secondary observables is shown in fig. 8. While an 
excess of events over both DIS Monte Carlos is observed also here, differences 
in the shape of the distributions are recognizable; neither the shape of the 
DIS Monte Carlos nor the shape of the signal expected from QCDINS are 
well reproduced by the data. 

These results do not yet allow any firm conclusion on a possible contribution 
of /-induced processes to deep inelastic electron proton scattering. Never- 
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Figure 7: Observables used to cut ("primary observables" ) after cuts. From Ref. 

[16] 

theless, the observed excess of events and the gross similarity between data 
and predictions are intriguing and motivate the continuation of these studies 
with higher statistics and further refined techniques. The ultimate signal of 
/-induced effects would consist, of course, in demonstrating the violation of 
helicity-conservation in DIS-events. 
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